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Abstract 

Let S{n, k) denote the Stirling number of the second kind, and let Kn be such that 
S{n,Kn-l) < Sin,Kn) > 5(n,if„ + l). 
Using a probabilistic argument, we show that, for all n > 2, 

Le"'(")j - 2<Kn< Le""^"^] + 1, 
where [x] denotes the integer part of x, and ■w{n) denotes Lambert's W function. 

1 Introduction 

The Stirling number of the second kind, denoted S{n,k), plays a fundamental role in many 
combinatorial problems. It counts the number of partitions of {l,...,n} into k non-empty, 
pairwise disjoint subsets, and may be defined recursively as 

S{n, k) = S{n-l,k-l) + kS{n -l,k), n>l,k> 1, 

together with 5(0,0) = 1, 5(n,0) =0, n > 1. 

According to Harper for each n > 1, the polynomial X^^^q S{n, k)x^ has only real zeros. 
By Newton's inequalities ([7], p. 52), log5(n. A;) is strictly concave in k. It follows that there 
exists some 1 < Kn < n such that 

5(n, !)<...< 5(n, Kn) > S{n, K„ + 1) > . . . > S{n, n). 



In other words, the sequence S{n,k), k = l,...,n, is unimodal, Kn being a unique mode if 
S{n,Kn)^S{n,Kn + l). 

Determining the value of Kn is an old problem ([SlIiniEllIlEllIllllSKIlllS]). A related long- 
standing conjecture ( [El O [12] ) is that there exists no n > 2 such that S{n, Kn) = S{n, Kn + 1). 
See [3j for a historical sketch and recent developments. 

In particular, Canfield and Pomerance [3j noted that 

K„G{[e»(«)j-l, [e"'(")j} (1) 

for both 2 < n < 1200 and n large enough (no specific bound is known on how large n has to be; 
see also [2]). Here and in what follows, [xj denotes the integer part of x and w(n) is Lambert's 
W function defined by 

n = u;(n)e"'('^). 

Based on this, it seems likely that ([1]) holds for all n. The purpose of this note is to present the 
following non- asymptotic bounds. 

Theorem 1. 

Le'^Wj - 2<Kn< [e'^^'^^J +1, n > 2. (2) 
Theorem [T] can be compared with the non-asymptotic bounds of Wegner |15j : 

Tl 

< 1 ^—^ ' ^ ^ 3; (3) 

log n — log log n 

n ( log log n — 1 \ 

Kn > 1 + , n > 31. 4 

logn \ logn / 

Note that the upper and lower bounds in ([2]) differ by 3, whereas the difference between the 
upper bound ([3|) and the lower bound ([4|) tends to oo as n — > oo. More precisely, it can be 
shown (details omitted) that the upper bound in ([2|) implies ([3|) if n > 7, and the lower bound 
in ^ imphes ^ if n > 34. 

In Section 2 we prove ([2]) using a probablistic result of Darroch [4| . The possibility of further 
refinements is discussed in Section 3. 
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2 Proof of dl) 

Recall Dobinski's formula 



°° un k 

e-Y,S{n,k)x^ = Y,^^ ^>1- (5) 



In particular 



Dividing ([5]) by ([6]) we get 



k\ 

k=l k=l 



e^S{n,k) = Y^-. (6) 

k=l k=l 



k\\^ ■^(^'^) fc _ 



-X' 



A: 



This has the following interpretation. If we let S be a random variable with probability mass 
function (pmf) Pr(S' = k) = S{n,k)/Y^'^=iS{n,i), k = l,...,n, and let Z he a Poisson(l) 
random variable independent of S, then the pmf of 5" + Z is 

Pr(5 + Z = fc)= ^ , k=l,2,... 

Ei=i ^ A' 

While the mode of S is hard to determine, that of 5" + Z is straightforward. (As usual, we 
call a random variable X on {0, 1, . . .} unimodal if its pmf is unimodal, and call any mode of 
the pmf a mode of X.) To relate the mode of S to that of S + Z, we invoke a classical result of 
Darroch |^ (see Pitman's survey [H]). Note that S can be written as a sum of n independent 
Bernoulli random variables since the polynomial X]fc=i ^{''^^ k)x^ has only real zeros. 

Theorem 2 ([J)- Let Xi, i = 1, . . . ,n, be independent Bernoulli random variables, i.e., each 
Xi takes values on {0, 1}. Then for any mode m of S = Y2^=i -^i 

\m-ES\ < 1. 

As a consequence of Theorem [2l we have 

Proposition 1. Let S = Y17=i -^i ^ '^^"^ '^f independent Bernoulli random variables. Let Z 
be a Poisson(l) random variable independent of S. Assume S + Z has a unique mode mi, and 
denote any mode of S by mo. Then 

fno < < iTT-o + 2. (7) 
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Proof. Note that, since the pmfs of S and Z are both log-concave, the pmf of 5* + Z is 
log-concave and hence unimodal. Denote /i = ES. By Darroch's rule, \fi — mo[ < 1. We show 
that Darroch's rule applies to S + Z, i.e., j/i + 1 — mi\ < 1. The claim then readily follows. Let 
Zk, k > 2, he Binomial(A;, 1/k) random variables, independent of S. Then S + Z^. is a sum of 
independent Bernoullis for which Darroch's rule applies; if we let be a mode of S + Z^, then 
— TTifcl < 1. Moreover, assuming mi is the unique mode of S+Z, we have lim^-^oo "m^k = fni- 
Thus 1^ + 1 — mi\ < 1. □ 

On the other hand, we have 

Proposition 2. For n >2, the sequence k"'/k\, k = 1,2, . . . , is unimodal with a unique mode 
at either k = [e'"^")] or k = [e'"^")] + 1. 

Proof. Denote u = e"'^"'^ and consider the ratio 

^ (fc + l)"/(fc + l)! ^ + 
^ k^/k\ k^ 

It is easy to see that f{k) ^ 1 for all k > 1. We also show that f{k) > 1 if k < u — 1 (i.e., 
k < \u\ — 1) and f{k) < 1 for k > u (i.e., k > \u\ + 1). The claim then follows. 

Noting that f{k) decreases in k, we only need to show f{u — 1) > 1 and f{u) < 1. However, 
direct calculation gives 

log /(n - 1) = -w{n) - n log (l - e""'^")) 
> -w{n) - n (^-e-"'(")) = 0; 
log f{u) = n log (l + e~"'(")) - log (e"'^") + l) 
< ne-'"^") - log (e"'(")) = 0. □ 

Then we obtain ([2]) as a consequence of Propositions [T] and O 

Corollary 1. Let n > 2, and denote K = [e^'^^^J. // k'^jkj. > {k^ + 1)"/(A:* + 1)!, then 
k-f — 2 < Kn < A;*; otherwise /c* — 1 < Kn < /c* + 1. At any rate ^ holds. 

3 Discussion 

A natural question is whether Corollary [1] can be further improved using this argument. This 
leads to an investigation of the bounds in ([7]). It turns out that the lower bound in ([7]) is 
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achievable. For example, in the setting of Proposition [H if we let n = 2 and Pr(Xj = 1) = 
1 — Pr(Xj = 0) = Pi, i = 1,2, with pi = 1/3 and p2 = 2/5, then mo = mi = 1 by direct 
calculation. It seems difficult, however, to find an example where the upper bound in ([7]) is 
achieved. After some experimentation we suspect that this upper bound is not achievable. This 
is further supported by the fact that, in the setting of Proposition [H we always have mi < mg + l 
when n < 5. To show this, let a = Pr(S' = i), i = 0, 1, . . .. By Newton's inequalities 



o (i + 2)(n-i) 

cj+i > r ■ n C,Q+2, < i < n - 2. 

[I + l)[n — I — 1) 



When n < 5 and 0<z<n — 2we have 

(i + 2)(n-i) 



> 2. 



{i + l)(n 

Thus cf_^i > 2ciCi+2 and c^q+i > 2cmoCmo+2 in particular. Since mo is a mode of S, we know 
Cmo > Cmo+i- Thus 

However, a simple calculation gives 



mo 



e[Pi{S + Z = mo + 1) - Pr(5 + Z = mo + 2)] = ^ 



Ck 



Cmo +2 



^ "2 Cmo+2 > U, 

which rules out mi = mo + 2 under the assumption that mi is the unique mode of 5 + Z. 
Conjecture 1. In the setting of PropositionUl mo < mi < mo + 1. 
It is clear that Conjectured] implies a sharper version of ([2]) 

Le«'(")j - l<Kn< [e'"^''^\ + 1; 
this is tantalizingly close to proving ([T|) for all n > 2. 
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